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Abstract 

The objective of this article is to highlight certain features of a number of coefficients that 
appear in models of phenomena of transport in anisotropic porous media, especially the coeffi- 
cient of dispersion, the 2nd rank tensor , and the dispersivity coefficient, the 4th rank tensor 
dijki, that appear in models of solute transport. Although we shall focus on the transport of 
mass of a dissolved chemical species in a fluid phase that occupies the void space, or part of it, 
the same discussion is also applicable to transport coefficients that appear in models that de- 
scribe the advective mass flux of a fluid and the diffusive transport of other extensive quantities, 
like heat. The case of coupled processes, e.g. the simultaneous transport of heat and mass of a 
chemical species, are also considered. The entire discussion will be at the macroscopic level, at 
which a porous medium domain is visualized as homogenized continuum. 

1 The coefficient of dispersion, 

We consider the transport by the mechanism of dispersion of a solute in a fluid that occupies the 
void space of a porous medium, or part of it. The coefficient of dispersion, Dij, appears in the 
Fickian-type expression for the dispersive flux of a solute (e.g. Bear, 1961) 

dc 

Ji = -Dij—, i,j = x,y,z, or i,j = 1,2,3, (1) 

where J« denotes the ith component of the solute flux vector J (= solute mass passing through a unit 
area of fluid in the porous medium cross-section, per unit time), and c is the solute's concentration 
(= mass of solute per unit volume of fluid). Furthermore, we use the term flux as an abbreviation 
for flux density. In ([!]), and everywhere else in this article, Einstein summation convention is 
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applicable, unless the sum symbol is used. The coefficient Djj is a 2nd rank tensor, relating the 
vector J to the vector Vc. 

Equation ([T]) is valid for the general case of an anisotropic porous medium, with the isotropic 
medium as a special case. The dispersion coefficient is characterized by: 

1. In thermodynamics, the rate of entropy production, S, is related to the thermodynamic driving 
force, X, and the thermodynamic flux, Y, called "conjugated flux and force" (De Groot and 
Mazur, 1962), by S = YiX^. Here, the flux of the solute, J is driven by — Vc, which acts as 
a "driving force." In this case, the rate of entropy production can be expressed by 



X 



x X 



, dc , 
{ ~dx~. 



^ dc dc 

£0 ' or D «ai&- £0 ' < 2 > 



in which, Y = %J and X = -^Vc. The (dimensional) parameter, x, depends on the 
considered transport phenomenon. Hence, the Dy— matrix is non-negative definite. 

2. The D«- matrix is symmetric, i.e., 

D i:j = Dji. (3) 

The above two statements are consequences of the fact that X and J are thermodynamically con- 
jugated force and flux (De Groot and Mazur, 1962), i.e., they satisfy 

8J L _dJ l 



dXj dXi ' 

We may comment here that the above considerations can also be applied to the hydraulic 
conductivity tensor, Ky, which is a 2nd rank tensor that appears in Darcy's law, for the general 
case of an anisotropic porous medium, say, 

Qi = (5) 
dx 5 

where qi denotes the ith component of the flux vector (= discharge per unit area of porous medium 
per unit time), h denotes the piezometric head, and the vector — V/i is the driving force. Here also, 
the Kjj— matrix is positive definite and symmetric, i.e., = Kjj. It is interesting to note that we 
may have cases with components K^- < for i ^ j. This is a consequence of the positiveness of the 
principal minor, 

K U K 22 - K\ 2 > 0, => -VKuK 22 < K 12 < y/KuKn. (6) 

For an isotropic porous medium, = K&y, with 6™ denoting the Kronecker delta. 

The above conclusions are valid also for the coefficient of thermal conductivity, the coefficient of 
(mass) diffusivity, and the dispersion coefficient of any extensive quantity transported in the fluid 
that occupies the void space, or part of it. For example: 
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1. Mass Transport of a dissolved solute. Such mass can be transported by advection 
with the moving fluid (not considered here), by dispersion (as described by ([I])) and by 
(molecular) diffusion. For a dilute system, the latter mass flux is described by Fick's law 
(for a porous medium), 

Jl ~ V dxj' [ ' 

in which denotes the ith component of mass flux of a dissolved 7-species (mass of 7 per 
unit area of porous medium per unit time), c denotes the concentration of that species, and 
T>1* denotes the coefficient of molecular diffusivity of 7 in the porous medium. It takes into 
account the coefficient of molecular diffusion of 7 in the fluid as a continuum, the porosity 
of the porous medium and the tortuosity of the void space within the latter (e.g. Bear and 
Bachmat, 1990, p. 193). The matrix Vj- is symmetric and positive definite. 

2. Heat transport. Although heat can be transported also by the solid matrix, we focus on the 
case in which the latter is thermally an insulator. Heat can then be transported by advection 
with the moving fluid (not discussed here), by thermal dispersion and by thermal conduction 
(which is a diffusive- type flux), described by (the averaged) Fourier law, 

dT 

in which jf 1 denotes the ith. component of the heat flux by conduction (heat per unit area 
of porous medium per unit time), T is the temperature, and A^* denotes the r/-component 
of the coefficient of thermal conductivity of the porous medium. Again, the latter takes into 
account the thermal conductivity of the fluid, as well as the porosity and tortuosity of the 
void space occupied by the fluid. The matrix A^* is symmetric and positive definite. 

When the fluid occupying the void space is moving, heat is also transported by advection, 
with the average velocity of the fluid, and by thermal dispersion. In analogy to (pQ), the 
thermal dispersive flux is expressed by 

rH* _ _p l H*^£_ (Q \ 

J di.S P ,i — U ij Q > W 

where J^ sp j denotes the ith component of the thermal dispersive flux vector (= heat per 
unit area of porous medium per unit time), and the vector — VT is the driving force. The 
D^*-matrix is symmetric and positive definite. 

Although, we have extended the conclusion about the coefficients being positive definite and 
symmetric, to a number of transport coefficients, there is a basic difference between the coefficients 



3 



of of hydraulic conductivity, K^, of mass diffusivity, 2?^- , and of thermal diffusivity in a porous 
medium, A|j*, as compared to those of (mass) dispersion, Dij, and thermal dispersion, Df-: 

• Kij depends on the geometry of the void space through which the water flows (tortuosity and 
width of pathways), in addition to fluid properties, such as density and dynamic viscosity. 

• T>J* and A^* depend on the geometry of the fluid pathways, as well as the diffusivity and 
thermal conductivity in the fluid. 

In all these cases, the geometry of the fluid-occupied domain is expressed by a scalar that 
represents the width of the pathways, and a 2nd rank tensor called tortuosity that represents 
the effect of the tortuous fluid pathways. Bear (1972, p. Ill) showed that the tortuosity is a 
2nd rank symmetric tensor. 

• The coefficients of mass and thermal dispersion, D,j and D^*, respectively, are functions not 
only of the geometry of the void space, but also of the velocity field within the porous medium 
domain. 

We have mentioned several times that a considered fluid may occupy only part of the void space, 
at some fluid saturation (= volume of fluid per unit volume of void space). Wherever a coefficient 
depends on the configuration of the void space, it also depends on the saturation of the considered 
fluid. This aspect is discussed in Sec. [H 

The last two coefficients are discussed in the next section. 



2 The dispersivity, a^ki 

The dispersivity, aijki, is related to the coefficient of dispersion Djj by (e.g. Bear, 1972, p. 610) 

U ij = a ijkl—^—, i-LUJ 

where V/% denotes the kth component of the fluid's velocity vector, V, with V = |V[. The coefficient 
a^jki is a fourth rank tensor having the following properties: 

1. Prom the expression for the rate of entropy production, S, and following the discussion leading 
to ([2]), we have 

dc \ ( dc \ 2 dc dc V^Vi 



S = Y,X, = x (- Dj ,_ j x * j = x%jU --— > 0, (11) 

from which it follows that is positive definite. Thus, all principal minors of are 
positive. 
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2. It has 3 4 = 81 components in a 3-dimensional porous medium domain, constrained by 2^ 2 " > 
1 = 511 inequalities. 



3. It is invariant under the following permutation of indices (see (1111) ) 

O'ijkl O'ijlki O'ijkl "•jikl- 



(12) 



Hence, only 36 components are independent of each other. Furthermore, there are 2 6 — 1 = 63 
constraining inequalities. 

Let us apply the above conclusions to an isotropic porous medium domain and to two anisotropic 
ones. 



2.1 Isotropic porous medium 

For this case, the 36 independent components reduce to two. This can be shown (Sirotine and 
Chaskolskaya, 1984, p. 651-2) by noting that the fourth rank tensor satisfies the relationships (|12p 
and is invariant under the the action of full rotational (orthogonal) symmetry group 0(3). In the 
case considered here, the Oy-jy-tensor can be represented by the matrix p a p, with 

a,/? =1,2,3 representing xx,yy,zz, a, (3 = 4, 5,6 representing yz,zx,xy, (13) 
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= o,t, we obtain: 



Pll - P12 



a-ijki = aT&ijdki H — ^—z — — (o~ik5ji + o~u5jk) 



and 



Dim = a T V5 im + (a L - a T 



VjVn. 

V 



(14) 



(15) 



(16) 



where the two independent coefficients are the longitudinal dispersivity, ai, and the transversal 
dispersivity, ax (e.g. Bear, 1972, p. 611). Furthermore, from the positive definiteness of a^i, it 
follows 

a L > 0, a T > 0. (17) 
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2.2 Axially symmetric porous medium 



This case, also called axially symmetric anisotropy, occurs, for example, when the porous medium 
is made of many, relatively thin isotropic porous layers. The same kind of material is obtained 
when we fill space by parallelepiped solid bodies, say boxes, a x a x b, in the x, y, z, directions, 
respectively, with equal spacing between all boxes. The z-axis is an axis of symmetry. In this case, 
the 36 independent components of aijki reduce to six (Sirotine and Chaskolskaya, 1984, p. 652), 

' Pll Pl2 P13 \ 
P12 Pll P13 



P31 P31 P33 
p 4 4 



P66 



Pll ~P12 



(18) 



p 4 4 o 
V Pm ) 

where we have made use of (|13p . taking the z-axis as the axis of symmetry. The number 6 is reached 
by considering fourth rank tensors that satisfy the relationships (|12|) and are invariant under the 
the action of the uniaxial symmetry group, D^h- 

Furthermore, because of the positiveness of the minors of the above matrix, the 6 constraints 
satisfy: 



' Pll Pl2 P13 ^ 



Pll,P33,PU > 0, Pll>Pl2, P11P33 > P13P31, det 



>0, 



(19) 



P12 Pll P13 
\ P31 P31 P33 J 

where we have made use of (113f> . 

Another form of presenting a^M, with the vector e (components ei) denoting the direction of 
the axis of symmetry, is 



a-ijki = aiSij5ki + -^-{SikSji + 5u5jk} + a^eiejSki + a 4 e k ei5ij + 
y {eie k 5ji + eje k 5u + eiei5 jk + ejei5 ik } + a§eieje k ei . 



(20) 



Poreh (1965), using Robertson (1940), derived an expression for a^ k i which missed 3 of the terms 
appearing in ([20]) . The corresponding Dij is 



DijV 



a\dij (V k V k )+a 2 V i V j + (V fe V fe ) + a 4 (V fc efc)(V r / eO-5 Ji + 
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-^{ei{V k e k )Vj + ej(V k e k )Vi + e^ViefiVj + ej(Viei)Vi} + a 6 eiej(y k e k )(Viei) 
eiej{a 6 {V k e k ) 2 + a 3 V 2 } + 6 tj {a A {V k e k f + ai V 2 } + a 2 V i V j + 
+a 5 (V k e k ){eiVj + ejVi}. 



(21) 



As an example, consider the case of a layered porous medium domain. This is an axially 
symmetric case, with the z-axis as axis of symmetry. Consider now uniform flow normal to the 
layers. This means: e\ = e 2 = 0, e% = 1 and V\ = V2 = 0, V3 = V. For this case, (l2T|) can be 
written as: 



V 



V, 



3THV 
*LVV 



a\ + a 4 , 

a\ + a 2 + 03 + a 4 + 2a 5 + a 6 



arw 
a T HV 
a L vv / 

For uniform flow parallel to the layers, say, in the +x-axis, e% = e 2 = 0, e% = 1 and V\ = V, V2 
V3 = 0, (j2Tj) can be rewritten as: 



(22) 



D 



V 



V, 



3LHH 
3THH 
3TVH 



CL\ + 02, 
Ol, 

a\ + a 3 . 



(23) 



3LH^ 
3 T HH 
3TVH ) 

Thus, to fully describe dispersion in a layered porous medium, when the flow is uniform and 
normal to the layers, we need one longitudinal and one transversal dispersivities. For uniform flow 
in the layers, we need (another) one longitudinal and two transversal dispersivities. In addition to 
these 5 coefficients, we need information on the direction of the axis of symmetry, i.e., we need six 
independent coefficients to completely describe the dispersion in such a domain, when the flows are 
uniform, wither along the axis of symmetry or normal to it. In the case of general flow through such 
domain, as shown earlier, we need six independent coefficients, plus information about the axis of 
symmetry. Note that because of the axial symmetry, only one angle is required to completely define 
the axis of symmetry. 

It is interesting to note that Batchelor (1946), in his work on axisymmetric turbulence, using a 
method based on invariance, introduced by Robertson (1940), derived a general expression for Dij 
that is based on five tensorial terms, e^ej, 5ij, ViVj, V^ej, Vjet, and four scalar functions, Cf. 



D i:j V = Ci(V, e)e iej + C 2 (V, e)<% + C 3 (V, e)V^ + C 4 (V, e){V t ej + V}eJ. 



(24) 



To construct the generic forms of the four scalar function Cj, where Cj = Cj(V, e),i = 1, ..4, 
such that they survive under the action of axisymmetric group, noting that there are only two 
quadratic-in-V scalar invariants: V 2 and (V^e^) 2 , we have: 



C\ = a 6 (y k e k ) 2 + a 3 V 2 , C 2 = a 4 (V fc e fc ) 2 + aiV 2 , C 3 = a 2 , C 4 = a 5 (V k e k ). 



(25) 



We have chosen the use of the same notation as in (I2ip in order to emphasize the explicit equivalence 
between the two expressions for Dij. Poreh (1965) followed the same analysis as that presented by 



Batchelor (1946). Altogether, the above comments explain why in the literature one can find the 
numbers four, five and six proposed by different authors. As we have shown, there is no contradic- 
tion; Djj is expressed by five tensorial terms, associated with four quadratic-in-F functions. The 
latter are constructed by six independent coefficients. 

Bear and Bachmat (1990, p. 215) also used Robertson's (1940) method and reached the results of 
2 independent coefficients for the isotropic case, and 6 for the axially symmetric one. However, they, 
erroneously, also assumed symmetry permutations of the couples of indices i, k and j, I, which led to 
only 5 independent scalar coefficients for the latter case. The reason for the additional assumption 
of symmetry followed from the assumption of analogy between dispersion and elasticity. 

2.3 Analogy with elasticity 

It is often stated that the relationship (fl~0j) is analogous to the relationship between the stress ((Jij) 
and the strain (e^) for an elastic solid: 

Sij = Cijj-i °~kh (26) 

in which Cijki denotes components of the solid's elastic modulus, which is a 4th rank tensor. Here, 
following the discussion leading to ([2]), the rate of entropy production, is expressed by: 

S = (xvij)(X£ij) = X 2 o-m {CijkiVij) , (27) 

which leads to the symmetry = C^uj, in addition to the validity of (|12|) . Thus, for an isotropic 
solid, have 

C^ki = Ai(SikSji + SuSjk) + A 2 5ij5ku (28) 

in which A\ and A 2 are scalar coefficients, which are related to the Lame constants for the elastic 
solid. Thus, for an isotropic solid, because these two coefficients are properties of the solid material 
only, they can be expressed in terms of the two coefficients: Young's modulus of elasticity and 
Poisson's ratio. 

In the axially symmetric case (symmetry group D^h), the total number (= 21) of independent 
components, Cyfc/ is reduced to 5. However, in the case of the dispersivity tensor, a^-jy, we have 
shown that the number of independent coefficients is 6 and not 5. Thus, the assumption of analogy 
is erroneous. 

2.4 Anisotropic porous medium with tetragonal symmetry 

We consider the case of an anisotropic porous medium with tetragonal symmetry. This kind of 
symmetry is referred to as (Sirotine and Chaskolskaya, 1984, p. 651). As an example of such 
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(29) 



porous medium material, we may consider one that is made up of orderly packed solid boxes ax axe 
with equal spacing between the boxes in all directions (or cubes with 3 different spaces). For this 
case, the dispersivity tensor can be described by the matrix representation 

/ Pll Pl2 Pl3 \ 

P12 Pll P13 

P31 P31 P33 

p 44 

p 55 
V p 55 J 

where we have made use of (|13j) . 

In this case, the 36 independent components can be expressed by 7 independent parameters. 
Furthermore, because of the positiveness of the minors of the above matrix, we have the 10 con- 
straints: 



' Pll P12 P13 ^ 



PU,P33,PM,P55 > 0, Pll>Pl2, P11P33 > P13P31, det 



> 0. 



P12 Pll P13 
\ P31 P31 P33 ) 

This is an obvious extension of (|39p . which presents the case of flow parallel to the layers in an 
axially-symmetrical case. Altogether, we need information on the directions in which the boxes, 
ax axe, are positioned in space (and this requires information on 2 angles), and then 3 coefficients: 
(one longitudinal dispersivity coefficient and two different transversal ones) for uniform flow in the a- 
direction, and 2 coefficients (one (different) longitudinal dispersivity and one (different) transversal 
dispersivity) for uniform flow in the c-direction Altogether, 7 independent coefficients (or 5, if the 
spatial directions are known). 

Another (non-matrix) form of the dispersion tensor, Dij, can be obtained by considering the 
invariants that can be constructed of the four vectors: the velocity, V and the three mutually 
orthogonal unit vector, a, (3, 7, such that 



a-kPk = Pklk = IkOLk = 0, a k a k = PkPk = Iklk = 1, 

7 = a. x (3, (3 = 7 x a, a = /3 x 7. (30) 
We also require that the symmetric tensor Dy (= Dji) survives under the action of tetragonal group, 



Dy(a) = Dy(-a), Dytf) = 0^-0), D y ( 7 ) = D^-7), Dy(a, 0, 7) = D 4J (/3, a, 7). (31) 
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With the above in mind, we start by listing all symmetric 2nd rank tensors built of the four vectors: 
V and a, (3, 7: 

Sij, caaj, Pifij, 7i7j, ViVj, (32) 

aiVj + ajVi, fcVj+foVi, jiVj + jjVi. (33) 

atiPj + ctjPi, 0aj + (3jji, ^aj+jjai. (34) 

Actually, because of the relationship (f30l) . only two of these vectors suffice to completely define 
D{j, provided, these vectors are chosen consistently. Here, once we select cx,f3, the vector 7 is well 
defined. Accordingly, we obtain 

S^, otiOLj, 0i(3j, ViVj, (35) 

aiVj + atjVi, (3iVj + /3jVi, cti/3j + ay/?j. (36) 

Next we have to select coefficients such that when multiplying each of the above terms, con- 
tributes a quadratic in V term to the expression of Djj. Furthermore, the resulting Djj has to 
satisfy (13ip . With this last requirement, the terms in (134j) have no contribution. The first three 
terms in (|35p . require coefficients of the form 

Ai {(V k a k ) 2 + (V k (3 k f) + A 2 V 2 , 

in which the ^4j's take on different values for each of the first three terms in f)35[) . For the last term 
in (|35p . the coefficient is a constant. As for the last term in (136p . it can be shown that only one 
scalar invariant can serve as a coefficient, as otherwise, condition (|3ip will be violated. Thus, we 
obtain: 

B {(a k V k ){aiVj + aj V^ + ((3 k V k )((3iVj + faVtf) . 
As to the last term in f|36j) . its contribution is: 

C(a k V k )(f3 k V k )( ai (3 3 +ajPi). 

Altogether, we arrive at the final expression, in terms of 7 coefficients: 

VDn = {A 1 {{a k V k f + {l3 k V k ) 2 )) + A 2 V 2 ) 5 i3 

+ {A 3 {(a k V k ) 2 + {f3 k V k ) 2 )) + A 4 V 2 ) (a iaj + (3^) 
+A 5 ({a k V k )( ai Vj + aj Vi) + (l3 k V k ){^V 3 + (3^)) 
+A 6 (a k V k )(f3 k V k )(a i (3 j + a.AO + A 7 V t V 3 . (37) 
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As an example, consider the case of a porous medium made of boxes, a x a x c, with a in the x 
and y directions and c in the z-direction, and a = (1, 0, 0), (3 = (0, 1, 0), 7 = (0, 0, 1). For uniform 
flow in the x or y directions, i.e., V\ = V, Vi = V3 = 0, (f37l) can be written as: 



A 



CLTHH 




\ 



arvH J 



V. 



O-LHH 
a THH 

arvH 



Ai + A 2 + A 3 + A 4 + 2A 5 + A 7 , 
Ai + A 2 + A 3 + A 4 , 
Ai + A 2 . 



(38) 



ij 



\ 



a THV 
O.LVV 



A 2 + A 4 , 
A 2 + A 7 . 



(39) 



For uniform flow in the +z-axis, V\ = V2 = 0, V3 = V, f)3T|) can be rewritten as: 

( aTHV ^ 
arnv V, 
a L vv J 

Thus, to fully describe dispersion in a layered porous medium, when uniform flow takes place 
in the +x or +y directions, we need one longitudinal and two transversal dispersivities-in the 
horizontal direction and in the vertical one. For uniform flow in the vertical (+z-)direction, we need 
(another) one longitudinal and one transversal dispersivities. Altogether we need 5 coefficients, plus 
information — two angles — on the directions of the axes, i.e., we need seven independent coefficients 
to completely describe the dispersion (i.e., Dy) in such a porous medium. 

Needless to add that in all cases, the actual values of the various coefficients have to be deter- 
mined experimentally. The knowledge about the number of independent coefficients in each case 
facilitates the design of the experiments. 

The analysis presented above for the cases of isotropic and anisotropic porous media, can be 
extended to any other types of anisotropy. 



3 Coupled processes 

3.1 Conjugate fluxes and forces 

The generalized Newton's law that describes molecular flux of linear momentum, Fourier's law, 
that describes conductive heat flux, and Fick's law that describes diffusive mass flux of a solute (= 
molecular diffusion) , are particular cases of the general linear law for phenomena of transport in a 
fluid phase (e.g. Bear and Bachmat, 1990, p. 99) 

Jf = -^L^^, i,j = 1,2,3, (or =x,y,z), (40) 
j=i 1 

where Jf denotes the zth component of the flux of an extensive quantity E n of a phase, $ n is a 
state variable associated with E n and L^J 1 is a coefficient of proportionality. 
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The three linear diffusive flux laws mentioned above, state that a nonuniform distribution of the 
state variable, <E> n (e.g., temperature), produces a flux of only the corresponding extensive quantity 
(e.g., heat). However, experimental evidence suggests that, in principle, gradients of state variables, 
<3? r , r ^ n, corresponding to other extensive quantities, may also contribute to the flux of E n . Such 
phenomena are referred to as coupled phenomena. An example of such phenomenon is the Soret (or 
thermodiffusion) effect, in which mass flux of a solute in a liquid phase is produced by a temperature 
gradient, in addition to the flux produced by the gradient of the solute's concentration according to 
Fick's law. Another example is the Dufour effect, in which heat flux is caused by a concentration 
gradient, in addition to the heat flux caused by temperature gradient, according to Fourier's law. 
Actually, the diffusive flux of a component, as expressed by Fick's law, manifests interdependence 
between the concentrations of components in a multicomponent system, 

a=l j=l 3 

where ./V denotes the number of dissolved chemical species, and c 13 denotes the concentration of 
the /3-species. Using the terminology introduced above, we have to regard the mass of a dissolved 
/3-species as one of the E n extensive quantities. 

All the above diffusive flux expressions are valid within the fluid that occupies the void space 
of a porous medium, or part of it, i.e., at the microscopic level of description. By appropriate 
averaging, we obtain the corresponding flux laws at the macroscopic, or averaged level. In fact, 
they have the same linear form as (|40p . except that J™ denotes the zth component of the flux 
of an extensive quantity E n of a phase, per unit area of porous medium, <3? n is the macroscopic 
state variable associated with E n , and the coefficient L^J 1 takes into account the porosity and the 
tortuosity of the porous medium. For short, we shall leave the same symbols also for the macroscopic 
interpretation, as our discussion focuses on forces, fluxes, and coefficients at the macroscopic level. 

We consider I extensive quantities, and t corresponding fluxes produced by the thermodynamic 
forces, 

d<$> r 

X] = - — , r = l,2,...,£, (42) 

Each diffusive flux Y[ and the corresponding force X\ are chosen such that they are thermodynam- 
ically conjugate force and flux (Sec. [I]). For the set of coupled fluxes, we have: 

e 3 

^ = EE^4 i = h2,3; q = l,2,...,t (43) 

r=l j=l 

From the phenomenological equations (|43p . it follows that the flux of an extensive quantity, E q , 
is a single- valued function of all the (coupled) thermodynamic forces associated with E q . They 
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express linear relationships between fluxes and thermodynamic forces. 

Equation (|43f) can be extended to coupled phenomena of a higher order, i.e., when the transport 
is described by coefficients that are tensors of 3rd and higher ranks. As examples, we may mention 
the coupling between piezoelectric and piezomagnetic phenomena discussed by Fel (2002). 

Before continuing with the discussion of coupled processes, let us make the following important 
comment. Although the hydraulic conductivity K%j, appearing in Darcy's law, and all the coeffi- 
cients L^j appearing in the equation that describes coupled diffusive fluxes, are 2nd rank tensor 
coefficients, there is a basic difference between them. Although Darcy (1856) obtained his law as 
an empirical one, it is actually a simplified version of the averaged momentum balance equation of 
a fluid. This law exists only at the macroscopic (porous medium) level, describing advective flux 
in a porous medium domain, while all other diffusive flux laws exist in a fluid continuum, i.e. at 
the microscopic level. When averaged, they become macroscopic diffusive flux laws (of heat, mass 
of a solute, etc.) in a porous medium domain. Coupling as discussed above exists only among 
diffusive- type fluxes. 

In 1851, Stokes postulated that in (f40j) . the coefficients L|J are symmetric with respect to the 
permutation of the coordinates i and j, i.e. 

Li] = L« (44) 

This postulate states that the transformation of a unit force along one axis into a conjugate 
flux along another axis remains unaltered when those axes are interchanged. In Sec. [U we have 
discussed the three coefficients Djj, and D^*- 

Bear and Bachmat (1990, p. 226) have extended the concept of coupled processes to the macro- 
scopic level in a porous medium domain. They suggested that the concept of coupling, remains 
valid also at the macroscopic level, provided the coefficients involved, L?J corresponding to a phase 
continuum, be multiplied by the tortuosity of the porous medium. 

Let us focus on the cross coefficients for q ^ r, which give the flux of E q caused by the force, X r , 
associated with the gradient of e r (= the density of E r ). Employing the principle of microscopic 
reversibility of processes, and methods of statistical mechanics, Onsager (1931) showed that for the 
linear equations (|43p . and provided a proper choice is made for the fluxes, Y 9 , and forces, X r , the 
phenomenological coefficients are also symmetric in r and q, i.e. 

L% = m 1 * r. (45) 

This follows from the symmetry of the composite matrix. For the example, the composite matrix 
for 1 = 2 takes the form: 



13 



Lj 



r 22 



(46) 



/ 



These, relationships are known as Onsager's, or Onsager-Casimir's, reciprocal relations. They ex- 
press a relationship between any pair of cross-phenomena (e.g., thermal diffusion and Dufour effect) 
arising from simultaneously occurring irreversible processes (e.g., heat conduction and molecular 
diffusion). According to Onsager, the reciprocal relations, (|4"4"|) . hold under two conditions (e.g., 
Bear and Bachmat, 1990, p. 103): 

1. The relationship between each individual flux and its conjugate thermodynamic force is linear. 

2. The fluxes, Y 9 , and their conjugate forces, X 9 , should be selected such that 

I 3 

q=l i=l 

where S is the rate of entropy production of the system (De Groot, 1963, Ch. 1). 

Equation (|47p implies that each Y q and its conjugate force, X 9 , must be of the same tensorial rank 
(not necessarily a vector, as indicated by the single subscript). 
By the 2nd law of thermodynamics 



(q,r,i,j) 



(48) 



Hence, a necessary condition for the validity of (|48|) is that all diagonal matrices L?J are positive 
definite. Let us consider three special cases of porous media. 



3.2 Isotropic porous medium, 0(3) 

For an isotropic medium, using the example of i,j,= 1,2,3 and t = 2, e.g., coupled fluxes of heat 
and solute, we can write the coefficients L?J in the matrix form, 



/ / a N 
a 
0a 



V 



/ c \ 
c 
0c 



/ c \ \ 

c 
0c 



/ d \ 
d 
d 



(49) 
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with a,c,d denoting coefficients, e.g., a = L\\, and d = L 22 . We have here 3 independent coeffi- 
cients, with the inequalities a, d > 0, and ad > c 2 . 

3.3 Anisotropy — with uniaxial symmetry, 

For this case, again, with i, j, = 1, 2, 3 and £ = 2, the matrix of coefficients, can be represented as 

//a00\ / 6 \ \ 



a 
d 



6 
0c 



V 



/ 6 o \ 
6 
0c 



( e ^ 
e 

V° // 



(50) 



/ 



Altogether, we have here 6 independent coefficients, which are governed by 6 inequalities: 

a, d,e,f> 0, ae > b 2 , df > c 2 . 



3.4 Anisotropy— with orthorhombic symmetry, D 2 h 



For the example of £ = 3, the matrix takes the form, 

/ / a \ / m n n \ 
6 
0c 



V 



/ 



m 
n 
0s 



/ r N ^ 

p 
V q J 



( m N 
n 
0s 



V 



( r ^ 

p 
V q ) 



1 d A 

e 

v° ° /; 

/ u o o N 

v 

10 



u 






















w 



9 
h 
i 



(51) 



Here we have 18 independent coefficients, satisfying the 18 inequalities: 

a,b,c,d,e,f,g,h,i > 0, 



ad > m 2 , be > n 2 , cf > s 2 , ag > r 2 , bh > p 2 , ci > q 2 , dg > u 2 , eh > v 2 , fi > w 2 . 
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For any £, the number of independent components is: 

for symmetry group 0(3) : -£(£ + 1). 

for symmetry group AxA : £{i + !)■ 

3 

for symmetry group D 2 h ■ -£(£ + 1)- 

It it is interesting to note that the material symmetries, D^, and Dq^, with horizontal honey- 
comb cross-section, which are close Ax>/i have also £(£ + 1) independent coefficients. 

4 Multiple fluid phases 

So far, the discussion has been related to phenomena of transport in saturated flow, i.e., in a 
fluid that fills up the entire pore space. However, often, a number of fluid phases occupy the 
void space simultaneously, each at at a certain saturation, where the saturation, S a of an a-phase 
is defined as the ratio between the volume of the a-phase and the volume of the pore space, 
with Yl a = 1- Under such conditions, considering phenomena of transport within a fluid phase, 
neglecting interphase transfers of the considered extensive quantities across (microscopic) interphase 
boundaries, each of the components: Kij, aijki, and all Ljps, is a function of the phase saturation, 
S a . The reason for this dependence stems from the fact that the geometric features of the phase 
occupied portion of the pore space, e.g. the tortuosity, varies with the saturation. 

For example, for an anisotropic porous medium, with orthorhombic symmetry and with x, y, z 
principal directions, k xx = f x (S w ), k yy = f y (S w ), k zz = f z (S w ), but still, the K xy = f xy (S w ) = 
K yx (S w ), etc. Another example is the dispersivity. For an isotropic porous medium, we have 

ah = faL(S w ), a L = faT(S w ), with f aL / f aT . 

5 Summary 

We have studied a number of transport coefficients that appear in the laws that govern advective, 

diffusive and dispersive fluxes of extensive quantities that are transported within a phase occupying 

the void space of a porous medium or part of it. We focused on an anisotropic porous medium. 

Making use of the basic features of tensors of 2nd and 4th rank, we have determined for each 

considered transport coefficient, the number of independent components. This information will 

enable further research on the relationships among the various dispersivity components. It will also 

facilitate the analysis of experiments that are conducted in order to determine the numerical values 

of these components by inverse methods. 

In a separate article, Bear and Fel (2009) discuss the minimum number of experiments that is 

required in order to determine the values of the six a« moduli for the axially symmetric case. 
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